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Abstract. Giant resonances in the continuum for deformed nuclei are studied with the time-dependent 
Hartree-Fock (TDHF) theory in real time and real space. The continuum effect is effectively taken into 
account by introducing a complex Absorbing Boundary Condition (ABC). 

PACS. 21.60.Jz Hartree-Fock and random-phase approximations - 24.30.Cz Giant resonances 
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Giant resonances (GR) are the most prominent fea- 
ture of nuclei at high excitation energy. Recently, the GR 
becomes an active issue again in relation with the struc- 
ture of unstable nuclei. The major part of GR's strength is 
embedded in the continuous energy spectrum, at energies 
higher than the neutron and proton separation energies. 
Thus, the coupling to the continuum is a key issue to un- 
derstand the excitation mechanism and damping of the 
GR. 

It has been well-known that the effects of the con- 
tinuum can be properly treated in the random-phase ap- 
proximation (RPA) with Green's function in the coordi- 
nate space [Q. However, it is very difficult to directly ap- 
ply the method to deformed nuclei because construction 
of the Green's function becomes a difficult task for the 
multi-dimensional space. Recently, we have developed a 
numerical method of constructing the three-dimensional 
(3D) Green's function with a proper boundary condition 
and studied photoabsorptions of molecules ||]. Although 
the method works successfully in calculations of response 
functions in electron systems, complexity of nuclear Hamil- 
tonian makes its application to nuclear physics harder 
1^. Therefore, we investigate an alternative treatment of 
the continuum in deformed systems; Absorbing Boundary 
Condition (ABC) method. The ABC method was exten- 
sively utilized in a field of the quantum chemistry, in order 
to study chemical reactions (for instance, see Ref. [|). We 
have shown that the ABC method is also very useful to 
treat the electronic continuum in molecules and clusters 
1^. We are also applying the ABC method to studies of the 
continuum effects in nuclear response [^,^ and in nuclear 
reaction 

The essential trick for the treatment of the contin- 
uum in the ABC method is to allow the infinitesimal 
imaginary part in the Green's function, ie, to be a func- 
tion of coordinate and finite, ie(r) Q. The e(r) should 



be zero in the interaction region and positive outside the 
physically relevant region of space. This is equivalent to 
adding the absorbing potential to the original Hamilto- 
nian, H H — ie{r). We adopt a function form of e(r) 
being linearly dependent on the radial coordinate r. The 
conditions and limitations on e(r) are discussed in a num- 
ber of works (see Ref. and references therein). 

In this work, we use ABC in the time-dependent Hartree- 
Fock calculations on a 3D coordinate grid. In the real-time 
calculations, the linear response is computed by apply- 
ing an impulsive external field to the Hartree-Fock (HF) 
ground state, then calculating the expectation values of 
some observables as a function of time, and then Fourier 
transforming to get the energy response |^,|[||]. Since all 
frequencies are contained in the initial perturbation, the 
entire energy response can be calculated with a single time 
evolution. 

We use the Skyrme e nerg y functional of Ref. ||l^ with 
the SGII parameter set [|ll[. For calculations of the HF 
ground state, the imaginary-time method of Ref. |12 is 
utilized. For the time evolution of the TDHF state, to 
which the time-odd components in the energy functional 
also contribute, we follow the prescription in Ref. ||l^. The 
model space is a sphere whose radius is 22 fm. The absorb- 
ing potential,— ie(r), is zero in a region of r < 10 fm, while 
it is non-zero at r > 10 fm. The TDHF single-particle 
wave functions are discretized on a rectangular mesh in a 
3D real space. In order to reduce number of mesh points 
outside the interaction region, we employ curvilinear co- 
ordinates which are mapped by a change of coordinates 
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to an adaptive mesh in Cartesian coordinates 1 14 
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details will be published in our forthcoming paper 

In this paper, we discuss an application of ABC to the 
isovector giant dipole resonance (GDR) in ^''Mg. The den- 
sity distribution of the ground state has a prolate shape. 
In Table |l|, calculated single-particle energies of neutrons 
and protons are listed. The deformation lifts the (2j -f- 1)- 
fold but not the Kramers degeneracy. The energies of the 
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highest occupied single-particle orbitals for neutrons and 
protons are —14.4 MeV and —9.7 MeV, respectively. Since 
the height of Coulomb barrier for protons is 4 ~ 5 MeV, 
the nucleonic continuum plays an effective role when the 
nucleus gets excitation energy larger than 14 15 MeV. 
We take an initial external field as 

Ve.t(t)-fcr|i(l-r,)e-^|<5W, (1) 

where k should be small enough to validate the linear re- 
sponse approximation. Of course, after the Fourier trans- 
form, the calculated El strength should not depend on the 
magnitude of k. In this calculation, we choose k = 0.001. 

Fig. (a) shows that the time evolution of El dipole 
moment initiated by the impulsive perturbation of Eq. (|I|) , 
{Vext)x- We perform two kinds of TDHF simulation with 
the same initial perturbation; TDHF with the ABC de- 
scribed above, and the one with the box (vanishing) bound- 
ary condition at r = 10 fm. Time evolution is carried out 
up to T = 10 ?i/MeV. The total period of time evolution, 
T, is related with the energy resolution after the Fourier 
transform, AE, by AE — 2nT%/T. In the box calculation, 
the dipole oscillation damps in a short period of time. 
However, we find a beating structure in a long period. 
The beat comes from reflection of outgoing waves dur- 
ing the time evolution of the TDHF state. On the other 
hand, the ABC calculation does not indicate any promi- 
nent beating pattern. Since the difference becomes signif- 
icant aX t > 1 ?i/MeV, the box calculation can provide a 
reliable result only when the strength function is averaged 
over an energy width of 6 MeV. 

The Fourier transform of the dipole moment, Dx{t), 
Dy{t), and Dz{t), calculated with ABC leads to El os- 
cillator strengths of ^''Mg shown in Fig. |l| (b). Here, we 
use a smoothing parameter of = 1 MeV. The calcula- 
tion shows a prominent double-peak structure which cor- 
responds to the deformation splitting of the GDR. The 
distribution of the oscillator strength qualitatively agrees 
with the photoneutron cross section [Q. A high energy 
tail of the oscillator strength (at E > 25 MeV) cannot be 
reproduced in the calculation with the box boundary con- 
dition. It turns out that the box calculation significantly 
changes the El strength at £' > 23 MeV, compared to the 
ABC results of Fig. |] (b). In other words, the damping ef- 
fect caused by fast emissions of a nucleon is important for 
E>23 MeV. 

In summary, the isovector GDR in the continuum of 
deformed nuclei is investigated with the linear response 
theory. The single-nucleon continuum is taken into ac- 
count by the ABC. The real-time TDHF is an efficient 
method to calculate the energy response over a wide range 
of energy. We consider that the present method, TDHF 
combined with ABC, is promising and we are studying ex- 
cited states and responses in nuclei near drip lines P,|l^ . 
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Table 1. Occupied single-particle energies for ^*Mg calculated 
with SGII. Each level has a two-fold degeneracy. 
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Fig. 1. (a) Calculated El moment as a function of time, Dx(t), 
for ^''Mg with ABC (thick line) and with BBC (Box boundary 
condition: thin line). The direction of the initial external dipole 
field is taken to be perpendicular to the symmetry axis (z-axis). 
(b) El oscillator strength calculated by Fourier transforming 
D(t). 
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